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Like many other intermetallic materials, MoSi 2 coatings are typically hard, but prone to catastrophic failure 
due to their low toughness at ambient temperature. In this paper, a self-assembled hierarchical structure that 
closely resembles that of nacre (i.e., mother of pearl) was developed in a MoSi 2 -based coating through a 
simple, yet cost-effective, depostion technique. The newly formed coating is tough and can withstand multiple 
indentations at high loads. Key design features responsible for this remarkable outcome were identified. They 
include a functionally graded multilayer featuring elastic modulus oscillation, varying sublayer thickness and 
a columnar structure that are able to attenuate stress concentrations; interlocking boundaries between 
adjacent sublayers that improve the bonding and arrest the cracks; a transitional layer that bridges the coating 
and substrate and facilitates load transfer. Moreover, the contributions of six important structural 
characteristics to damage resistance are quantified using finite elemnet analysis and in an additive manner 
(i.e., from low- to high-level complexity). The in-situ toughened coating is envisaged to enhance the 
mechanical performance and extend the lifespan of metal components used in safety-critical applications. 



Applying coatings for surface protection can be traced back more than 20,000 years, when natural materials 
such as blood, eggs or plant fibres were used 1 . Despite such a long history, major technological advances 
only took place about half a century ago, when physical vapour deposition was utilized to prepare hard, 
adherent coatings, significantly enhancing the performance and durability of cutting tools, engine components 
and biomedical devices 2 . 

With the introduction of extremely hard coatings in recent years, catastrophic failure, resulting from their low 
tolerance to structural defects and damage, has become a limiting factor in safety- critical applications 3,4 . In 
contrast, a combination of high hardness and good toughness can be found in many naturally occurring materials, 
for instance, nacre, enamel and bone 5 " 7 . These materials exhibit a multilayered hierarchical structure, built from 
mineral crystallites and soft organic compounds 8 . Multiple toughening processes work together to make these 
biological materials exceptionally resistant to contact damage 9 . In addition, biological interfaces, such as enamel - 
dentin junctions, exhibit interlocking traits 10 , offering distinct advantages in facilitating load transfer and enhan- 
cing interfacial strength 11 . 

Over the past several years, considerable efforts have been made to develop hard, yet tough, materials and 
coatings, often through emulating materials design in nature and backed by theoretical analysis and simula- 
tions 12 " 17 . Multilayered coatings modelled on nacre that comprise alternating hard and soft layers were reported 
and attracted considerable attention 18 " 21 . They were commonly fabricated using an additive process, i.e., sequen- 
tial deposition of hard and soft materials under precisely controlled deposition conditions. The sophisticated 
procedures required for synthesising these coatings render manufacturing processing prohibitively expensive for 
many industrial applications. Additionally, the additive fabrication makes the interfaces potentially weak and 
vulnerable to fatigue damage. Moreover, unlike their natural counterparts, multiple toughening processes are 
rarely seen in these synthetic coatings. 

Recently, we have developed a bimodally grained microstructure in a MoSi 2 -based coating by using a simple, yet 
cost-effective, fabrication technique - glow discharge plasma deposition 22 . The resultant coating displays an 
increased resistance to contact damage through crack bridging and deflection. Driven by this advance, a complex 
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hierarchical structure was created in situ for a MoSi 2 -based composite 
coating in this study. This newly developed coating captures key 
design elements of tough biological materials and, as such, possesses 
a remarkable damage tolerance as well as excellent load- carrying abil- 
ity. Important structural factors responsible for the improved tough- 
ness were modelled and evaluated. The work reported here opens a 
new pathway to developing hard, yet tough, coatings for safety-critical 
applications, where catastrophic failure cannot be tolerated. 



Results 

Unlike nanocrystalline structures reported for monolithic MoSi 2 and 
Mo 5 Si 3 coatings 22 ' 23 , in this study a complex multilayer structure was 



developed in-situ for a MoSi 2 -based composite coating (named here 
£ MMC). To reveal its microstructure features, both plan-view and 
cross-sectional images of the MMC coating were acquired. As shown 
in Fig. la, the coating contains three distinct layers. The outermost 
layer is — 1 um thick (refer to the top surface region labelled as 
"MoSi 2 " in Fig. la) and comprises fine MoSi 2 crystals with an average 
diameter of ~5 nm, according to plan- view TEM shown in Fig. lb. 
The fringe spacing of the crystallites, outlined by dotted red circles in 
Fig. lb, is 2.18 A, corresponding to the (111) lattice plane of hexa- 
gonal C40- structured MoSi 2 . These fine C40-MoSi 2 grains are 
arranged with a (111) texture (see XRD data shown in Fig. lb). 
Note, the Mo 5 Si 3 peaks in the XRD pattern are supposed to originate 
from region underneath the outer layer. 



• C40 MoSi. 




Figure 1 | (a) Cross-sectional bright-field transmission electron micrograph (TEM) revealing the pristine microstructure of the MMC coating, (b) X-ray 
diffraction pattern showing that MoSi 2 is the dominant phase within the outer layer of the coating. Note that the peaks for Mo 5 Si 3 originate 
presumably from region beneath the outer layer. Below the XRD curve is a plan-view bright-field TEM image taken —500 nm deep from the surface in 
region outlined by solid white line in a). (c,d) Enlarged field emission scanning electron microscope (FESEM) images of regions "c" and "d" outlined by 
dotted blue lines in a), (e) Cross-sectional TEM image of columnar structure of region "e" outlined by solid yellow line in d). (f) 2D elemental mapping of 
region "f ' outlined by dotted red line in a) . (g) Cross-sectional TEM image of the junction between the intermediate layer and substrate, as outlined by the 
dotted blue line "g" in a). 
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Figure 2 | (a) Schematic illustration of the self- assembled multilayer structure in the MMC coating (not drawn to scale), (b) Plan-view bright-field TEM 
image taken from the area —5 (am below the surface as indicated by dashed white horizontal line in a), (c) Plan-view bright-field TEM image taken from 
the area —10 (am below the surface as indicated by dashed yellow horizontal line in a) , with corresponding SAED pattern of "coarser" grains shown in the 
lower left inset. Note in a), below the outer layer ( — 1 (am thick MoSi 2 > refer to Fig. la) is the intermediate layer (—8.5 (am thick, as seen in Fig. la) 
exhibiting a graded, multilayer structure that comprises —10 periods (aslo refer to Fig. la). Each period consists of 2 sublayers; namely, the Mo 5 Si 3 layers 
(in blue colour, designated as "Sublayers A" in Table 1) and the MoSi 2 -rich columnar layers (in red and yellow, designated as "Sublayers B" in Table 1) 
with Mo 5 Si 3 grains embeded (in blue). 



Below the outer MoSi 2 layer is a —8.5 (am thick intermediate layer. 
It exhibits a multilayered structure, in which sublayers interlock with 
each other. Fig. lc and Id are enlarged views of multilayer regions "c" 
and "d" outlined by dotted blue lines in Fig. la. The sublayers lighter 
in contrast, as further revealed in Fig. le, comprise tightly packed, 
vertically aligned rods, whose boundaries are highlighted by dotted 
white lines; the individual rods have a diameter of —100 nm and 
contain highly organized, elongated grains. The plan-view bright- 
field TEM images of the lighter contrast sublayers, taken from 
regions about 5 (am and 10 (am below the surface (as indicated in 
Fig. 2a), are shown in Fig. 2b and 2c, respectively. Both regions 
exhibit a bimodally grained microstructure: finer grains have a size 
similar to those in the outer layer and are indeed C40- structured 
MoSi 2 , whilst coarser grains, according to the associated selected area 
electron diffraction (SAED) pattern (Fig. 2c, inset), represent the 
D 8m Mo 5 Si 3 phase. The size of Mo 5 Si 3 grains increases from —50 
to —70 nm with increasing depth, whereas the diameter of MoSi 2 is 
smaller than 10 nm and barely changes throughout the intermediate 
layer. 

To determine the phase composition of the darker contrast sub- 
layers, EDS elemental mapping of region "f ' outlined by dotted red 
line in Fig. la is displayed in Fig. If. It reveals that these darker 
contrast sublayers are rich in Mo and Si. Moreover, the ratio of Mo 
to Si is found to be — 1 .7, implying that these sublayers are most likely 
the Mo 5 Si 3 phase. According to its concentration profile in Fig. If, 



Mo in the intermediate layer exhibits uphill diffusion behaviour. 
Uphill diffusion is a common feature of multi-component systems, 
occurring when the activity gradient of a component is greater and in 
the opposite direction to its concentration gradient 24 . Uphill dif- 
fusion of the Mo element observed here may be explained as a con- 
sequence of a negative chemical interaction between Ti and Mo. 
When the Mo element diffuses toward the coating/substrate inter- 
face, the presence of Ti decreases the activity of Mo to the extent that 
the Mo activity gradient in the intermediate layer is in the opposite 
direction to the Mo concentration gradient. The phenomenon of 
uphill diffusion of Mo was also observed in the Me vs. MoSi 2 dif- 
fusion couples where Me = W, Re, Nb or Ta 25 . The transitional layer 
between the intermediate and the substrate is about 1 |am thick 
(Fig. la, region "g"); the main constituent in this region changes 
from fine Ti 5 Si 3 (the upper half in Fig. lg) to relatively coarse P-Ti 
grains (the lower half in Fig. lg). Below the P-Ti rich layer is a region 
(—7.3 (am thick) made of a" and a' phases (Fig. la). The type of 
constituents present in this region seems to be closely associated with 
the Mo contents. The upper part, having a Mo content greater than 
6 wt%, is composed of acicular a" phase, while the lower part is 
dominated by platelet-like oc' phase. 

According to the preceding observation, a sophisticated hierarch- 
ical structure that closely resembles that of nacre is developed in the 
MMC coating, as illustrated in Fig. 2a. The basic building blocks are 
MoSi 2 and Mo 5 Si 3 . The elongated MoSi 2 grains form rods embedded 
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with Mo 5 Si 3 ; each rod is —100 nm in diameter. They are bundled 
together into layers assuming a columnar structure (as seen in red 
and yellow colours in Fig. 2a), through which the content of Mo 5 Si 3 
increases. These layers are locked together by Mo 5 Si 3 sublayers (as 
seen in blue in Fig. 2a), and thus a functionally graded, multilayer 
structure results. Unlike those reported in the literature 26 ' 27 , such a 
nacre-like structure is formed in situ in this work. Some striking 
hallmarks describing the spinodal decomposition of alloys, such as 
periodic microstructure and compositional fluctuation, are observed 
in the MMC coating 28 . The spinodal phase segregation during the 
coating growth is a thermodynamically driven and diffusion -con- 
trolled process. To make it occur, high deposition temperature is 
introduced in this work. However, to understand the spinodal 
decomposition in the MMC coating, detailed thermodynamical ana- 
lysis is needed, for example, the determination of the Gibbs free 
energy of the mixed MoSi 2 /Mo 5 Si 3 system, which is part of ongoing 
investigation. 

Typical load- displacement curves of instrumented indentation 
testing on the MoSi 2 , Mo 5 Si 3 and MMC coatings are shown in 
Fig. 3a. The hardness (H) and elastic modulus (E) of the coatings, 
calculated according to the Oliver- Pharr method 29 , are given as a 
function of the maximum applied load (P) in Fig. 3b. The indentation 
size effect can be seen with lower P leading to greater H and E values. 
It is probably associated with the strain gradient and/or cracking 
underneath indenter. The H values of the MMC coating are greater 
than those of monolithic MoSi 2 and Mo 5 Si 3 coatings. However, the E 
values of the MMC coating are comparable to that of MoSi 2 at lower 
loads (20 and 50 mN), because the outer layer (~1 um thick) of the 
MMC is made of MoSi 2 . With the increase of indentation load (or 
depth), the modulus of the MMC is greater than MoSi 2 . To appreci- 
ate how the as-prepared coatings would behave when subjected to 
indentation loads, five Vickers indentations were performed in a 
pattern having one surrounded by four others (Fig. 4a-d). The dis- 
tance between the central indent and any of surrounding ones is 
about twice the indentation diagonal length. For the monolithic 
MoSi 2 and Mo 5 Si 3 coatings, both edge ring and corner radial cracks 
are visible at a load of 300 g (Fig. 4a and 4b). In contrast, for the 
MMC coating, except for ring cracks around the indents, no radial 
cracks are seen at the indent corners, even at 1000 g (Fig. 4c and 4d). 
The SEM image in Fig. 4e shows that the annular rings also appear 
within the indentation site. Such a deformation mode has been 
observed when indenting hard coatings having a columnar-grained 
structure 30 , suggesting that the plastic deformation might be enabled 
by shear sliding in the coating. To obtain a deeper picture of the 
damage resistance of the MMC coating, the indents created by a 
Berkovich indenter at 500 mN were sectioned and examined first 
using an FIB (Fig. 4f) and then by TEM (Fig. 4g). No open cracks can 
be detected in the subsurface of the MMC coating and the occurrence 
and arrest of shear sliding is evident ~3 urn beneath the indenter. In 
comparison, subsurface cracks were found to occur and propagate in 
the MoSi 2 and Mo 5 Si 3 coatings with little resistance (Fig. 4h&i). 
From these observations, the reason that the H and E of the MMC 
coating are greater than the MoSi 2 is probably due to the formation of 
cracks in the brittle MoSi 2 coating (refer to Fig. 4h) that lowers the 
measurement values. In addition, the Mo 5 Si 3 layers in the MMC 
coating serve as physical barriers that enhance the rigidity (i.e., 
modulus) and hardness of the coating by resisting the sliding of 
MoSi 2 grains. 

In an earlier paper on TiSiN-based systems by Wo et al. 2 \ two 
factors; that is, the shear sliding and modulus oscillation, were uti- 
lised in the design of tough coatings. Moving beyond that, numerous 
structural and mechanical features are identified in this work. To 
quantify the contributions of key features including the composi- 
tion/property gradient and the changing thickness of hard sublayers 
to the remarkable toughness of the MMC coating, finite element 
analysis was used to evaluate the stress concentration induced by 



20 



0 



M05S13 


a) 


MoSi : 




-*- MMC 





















0 125 250 

Displacement hi nm 




0 100 200 300 400 500 
Load P mN 

Figure 3 | (a) Load-displacement curves of the as-deposited coatings 
obtained using a Berkovich indenter at a maximum load of 20 mN. (b) The 
variation of indentation hardness (H) and elastic modulus (E) of the 
coatings at different loads. Ten indentations were made at each load, with 
the results presented here representing average values. Scatter bars 
representing ± one standard deviation are included in the data. 

indentation loading in coatings. Seven models with increasing struc- 
tural complexity, designated as M0 to M6, were constructed: M0 
represents a monolithic coating made entirely of MoSi 2 ; Ml also 
signifies a monolithic coating made of MoSi 2 except that a trans- 
itional layer is inserted between the coating and substrate; M2 - M6 
are multilayer coatings with distinct features in layer thickness, com- 
position and deformation modes (Table 1 and Fig. 2). The mech- 
anical properties of structural components - MoSi 2 and Mo 5 Si 3 - used 
in these models are given in Fig. 3b. 

Cracks in hard coatings are normally initiated by either shear or 
tensile stress. Considering this, the shear stress, t, induced by 
indentation to a depth of 2.4 um in the coatings MO-6 are shown 
Fig. 5. Here the shear stress is defined as (o n — g 33 )/2, where o n and 
o 33 are the first and third principal stresses, respectively. In the case of 
axial symmetry, t is within the r-z plane, and can be calculated using 
the equation t = {[{o rr — g zz )/2] 2 + T rz 2 } m . Notably, the maximum 
shear stress appears just below the contact surface, and is concen- 
trated within the MoSi 2 rich sublayers (due to its higher modulus) in 
the multilayer coatings (M2-6). A considerable reduction in both the 
stress magnitude and extent results, when the coating structure 
changes from monolithic to multilayered. This is particularly true 
for M6, which bears a striking resemblance to the nacre-like coating 
prepared in this work. For M6, the material volume populated by 
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Figure 4 | Scanning electron micrograph (SEM) of Vickers indentation patterns created on (a) Mo 5 Si 3 and (b) MoSi 2 at 300 g and on the MMC coating 
at (c) 300 g and (d) 1000 g. (e) SEM images of the edge of indents in Fig. 4d. Cross-sectional views of the indents created by a Berkovich indenter 
at 500 mN on (f) MMC, (g) also on MMC, showing the occurrence and arrest of shear sliding in a region —3 um deep from the surface, on (h) MoSi 2 and 
(i) Mo 5 Si 3 coatings. 



high shear stress, as shown in red in Fig. 5, is reduced markedly, 
compared with M0. Such a drop also occurs at other indentation 
depths modelled in this work. Further, the observed reduction is 
not only limited to the shear stress. Our modelling results indicate 
a reduction of all stress components within the multilayered struc- 
ture (i.e., M6), including the tensile stress and von Mises stress, 
compared to those in a monolithic coating (i.e., M0). 

Fig. 6 provides a quantitative comparison of the material volumes 
in seven types of coatings that are populated by various shear stresse 
levels. Notably, a consistent reduction in the stress is achieved by 
gradually increasing the structural complexity (i.e., incoporating new 
features to the coating design one step at a time). Judging from M0 
through to M6, the material volumes populated by shear stress level 



equal to or greater than 17.2 GPa are 17.9, 14.6, 14.4, 11.6, 8.8, 4.3 
and 1.2 |im 3 , respectively (refer to data in assocaited Table). In this 
case, adding a transitional layer (formation of new intermetallic 
compounds induced by interdiffusion between the coating and sub- 
strate) would reduce the stressed volume by —18% (Ml), while sub- 
sequent uniform insertion of alternating softer Mo 5 Si 3 layers (i.e., 
sublayers A in Table 1) with a thickness of 0.2 |im (M2) only brings 
about a marginal improvement (i.e., —1.4% over Ml). This may be 
due to a less pronounced modulus oscilation in the multilayer 21,31 . 
Advancing from M2, M3 features a gradual increase in the thickness 
of the Mo 2 Si layers (i.e., the sublayers B in Table 1), which results in a 
reduction of stressed volume by —35% over M0. Comparing with 
M3, the elastic modulus for the sublayers B is lowered in M4 from 
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to (E MoSi2 + E Mo5Si3 )/2 by varying compostion, yielding a 
reduction in affected volume by —51% over M0. For M5, the material 
volume associated with a stress 17.2 GPa or higher is reduced by 
—76% over M0, which demonstrates that a gradual change in the 
elastic modulus of the sublayers B across thickness (i.e., from E MoSi2 
to E Mo5Si3 ) is more effective than simply fixing the modulus of the 
sublayers B to the average value (as seen in M4). Finally, in the case of 
M6, in which the inelastic deformation occurs through shear sliding 
of the columnar elements in sublayers B, the extent of stress concen- 
tration is reduced markedly by —93% to only one-fifteenth that in 
M0. Furthermore, finite element calculation was conducted to quant- 
ify the load- carrying performance of the MMC coating. The results 
show that for the stress in the MMC to reach the same level as that in 
the monolithic coating (i.e., M0), about a 70% increase in applied 
load is required. This means that the MMC coating not only provides 
remarkable damage resistance, but also possesses excellent load- car- 
rying ability. 

The current trend in the development of nacre-mimetic materials 
and coatings is to mimic 'mortar-and-brick' structure of nacre. The 
resultant microstructure typically consists of hard ceramic platelets 
bonded together by polymer 1819,32 ' 33 . While such a design can 
enhance toughness, it is not suitable for applications involving 
extreme conditions, such as high speed machining where temper- 
ature may rise quickly under heavy load. To combat this challenge, 
the polymer component can be replaced by stronger, heat-resistant 
metals or ceramics 15 " 17 . In view of this, the MoSi 2 -based coatings 



comprising alternating layers of hard, thicker MoSi 2 and relatively 
soft Mo 5 Si 3 layers are developed in the present work. While the effect 
of elastic modulus oscillation may be less pronounced in this multi- 
layer structure, additional design features are introduced to ensure 
remarkable damage resistance can be achieved. 

Key toughening mechanisms enabled by the self- assembled hier- 
archical structure in the MMC coating can be described as follows: 
under initial contact, the functionally graded structure, created prim- 
arily by a continuous increase in Mo 5 Si 3 volume through the coating 
thickness, promotes the load transfer and, in so doing, discourage the 
stress build-up in the coating and at the coating/substrate interface. 
With an increase in load, deformation takes place through the plastic 
deformation of Mo 5 Si 3 components and/or the shear sliding in the 
columnar sublayers. The sliding process, once activated, can effec- 
tively mitigate the stress concentration without compromising the 
structural integrity of coating. For TiSiN-based coatings studied by 
Wo et al 21 , the columnar layers consist of vertically aligned TiN 
grains. In comparison, in the newly formed coating the columnar 
layers exhibit a hierarchical structure that comprises tightly packed 
rods; each has a diameter of —100 nm and contains highly orga- 
nized, elongated MoSi 2 grains. Consequently, shear deformation 
could occur at both 'rod' and 'grain' levels under high loads, making 
the new coating even tougher. It is also worth noting that Mo 5 Si 3 
sublayers have two functions: a) they serve as a physical barrier that 
resists the shear deformation of columnar sublayers, simultaneously 
increasing the stiffness, hardness and damage tolerance of coating; b) 



r/GPa 




Figure 5 | The distribution of the shear stress, r, obtained by FEA in (a) M0, (b) Ml, (c) M2, (d) M3, (e) M4, (f) M5 and (g) M6 at an indentation depth of 
2.4 urn. 



SCIENTIFIC REPORTS | 4:4239 | DOI: 1 0.1 038/srep04239 



6 




17.2 16.8 15.6 

Shear stress rl GPa 



14.4 




£/GPa 
380 



290 



200 

Figure 7 | Model configuration used in FEA. (a) The geometry of the 
model, (b) details of locally refined mesh elements in the coatings, and (c) 
the modulus distribution in the MMC coating. 
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Figure 6 | Material volume populated by different levels of shear stress in 
various coating structures (M0-M6) at an indentation depth of 2.4 urn. 

Data displayed in Table are an example showing material volume 
populated by stress stresses equal to or greater than 17.2 GPa. 

they interlock with columnar sublayers, giving rise to extended inter- 
facial area with strong chemical bonds, and are particularly effective 
in arresting lateral cracks. 

The relative thickness of soft and hard materials in a multilayer 
structure has an appreciable effect on its damage resistance. 
Theoretical analysis 34 and subsequent experimental observation 35 
have shown that hard, thicker layers alternated by thin, softer layers 
is most effective in resisting the crack propagation. Distinctly differ- 
ent from TiSiN-based systems studied earlier by Wo et al. 21 that 
features a less favourable design, i.e., soft thicker layers (TiN) sepa- 
rated by thin harder ones (TiSiN), here the new MoSi 2 -based com- 
posite coating exhibits a preferred modulus oscillation profile - hard, 
thicker layers (MoSi 2 dominated) separated by thin, softer layers 
(Mo 5 Si 3 ). Moreover, with an increase in depth, the thickness of hard 
MoSi 2 dominated sublayers increases, further enhancing the crack 
resistance. Last, but not least, a transitional layer is developed 
between the coating and substrate from the inter-diffusion of Mo 
and Ti elements. Its presence not only provides strong bonding that 
unites the coating with the substrate, but also further lowers the stress 
concentration at the coating/substrate interface. 

Discussion 

In summary, a MoSi 2 -based coating having a layered, nacre-like 
structure was synthesized in situ by using a simple, yet economic, 
fabrication technique - glow discharge plasma deposition. Inden- 
tation experiments show that the newly- developed coating possesses 
a remarkable resistance to contact damage. Finite element analysis 
demonstrates that effective lines of defence against damage exist in 
the coating, stemming from its self-toughened architecture. That is, 
a) a functionally graded multilayer, with elastic modulus oscillation 
that lowers the stress concentration that otherwise drives crack 
growth in the coating; b) the columnar substructure that enables 
inelastic deformation through shear sliding that dissipates energy; 



c) relatively soft Mo 5 Si 3 layers that enhance the bonding through 
interlocking and act as physical barriers against the shearing of 
columnar layers; d) the transitional layer that welds the coating 
and substrate together and promotes the load transfer through the 
coating/substrate interface. This work represents a notable technolo- 
gical advance in the synthesis and characterisation of bio -inspired, 
damage-resistant coatings that have potential to enhance the per- 
formance and durability of metal components operating under 
severe loading conditions. 

Methods 

Substrate discs, 40 mm in diameter and 3 mm in thickness, were cut from a com- 
mercially pure titanium rod (0.003% N, 0.010% C, and 0.074% O). Prior to coating 
deposition, the substrates were ground successively using SiC grinding papers of 400, 
800, 1200 grades and, finally, polished with 1.5 um diamond paste. The polished 
substrates were then ultrasonically cleaned in ethyl alcohol and dried in cold air. The 
MoSi 2 -based composite coating (termed here MMC) was deposited onto the sub- 
strate using a double cathode glow discharge apparatus, wherein a target having a 
stoichiometric ratio of Mo 3 oSi 7 o was mounted. For comparison, monolithic MoSi 2 
and Mo 5 Si 3 coatings were also deposited from targets having stoichiometric ratios of 
Mo 2 5Si 75 and Mo 5 oSi 5 o, respectively. The targets were fabricated from ball-milled Mo 
(99.99% purity) and Si powders (99.99% purity) by cold compaction under a pressure 
of 600 MPa. During deposition, one cathode was the target, and the other was the 
substrate 36 . The glow discharge sputtering conditions can be described as follows: the 
base pressure, 4 X 10" 4 Pa; target electrode bias voltage, —950 V; substrate electrode 
bias voltage, —400 V; working pressure, 35 Pa; target- substrate distance, 10 mm; and 
treatment time 3 h. The substrate temperature was set at ~800°C for the deposition 
of monolithic MoSi 2 and Mo 5 Si 3 coatings and at ~900°C for the MMC coating. 

The phase composition of the as -grown coatings was characterized with an X-ray 
diffractometer (XRD, D8 Advance, BRUKER AXS, Inc., Madison, WI 53711-5373, 
USA) operating at 35 kV and 40 mA. X-ray data were collected using a 0.1° step scan 
with a count time of 1 s. To reveal the microstructural features of coatings, plan-view 
samples for transmission electron microscope (TEM) observation were prepared 
using a single-jet polishing technique from the untreated side of the substrate 22 . The 
cross-sectional TEM samples were also prepared using a focused ion beam (FIB) 
microscope (FEI 200xP, FEI instrument, Hillsboro, USA) by sectioning through the 
area of interest. Detailed descriptions of the procedure can be found elsewhere 37 ' 38 . 
These samples were examined using a field emission gun TEM (Philips CM200, 
Eindhoven, Netherlands) operating at 200 kV. An energy dispersive X-ray spectro- 
meter (EDS) was interfaced to the TEM for elemental analysis. The crystal structure of 
the coatings was determined from selected-area electron diffraction (SAED) analysis. 
To investigate the damage patterns in the coatings, the surface and subsurface 
examination of indentations was performed using the FIB microscope and TEM. 

Instrumented indentation tests were conducted on the as-prepared coatings using 
a Berkovich indenter (NHT, CSEM Instruments, Switzerland). Fused silica was used 
as standard material in calibration. The load- displacement curves were obtained by 
driving the indenter with a constant loading rate of 40 mN/min into the samples at 
maximum loads ranging from 20 to 500 mN. The hardness and the elastic modulus of 
the coatings were calculated using the Oliver- Pharr method 29 . The damage tolerance 
of the as-deposited coatings was evaluated by a Vickers indenter over a range of loads 
up to 1000 g. 

Finite element analysis (FEA) was performed using COMSOL software to elucidate 
the roles of key microstructural features in mitigating the stress level and maintaining 
the structural integrity of the coatings. To do that a two-dimensional axisymmetric 
model of 100 X 100 um was constructed (Fig. 7a). It consists of a coating and 
substrate assembly that was loaded along the z axis by a diamond spherical indenter 
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with a radius of 10 urn. The thickness of the coatings is 10.5 |im. For coatings having 
a layered structure that mimics the MMC coating, they consist of an outer layer of 
1 um in thickness and an intermediate layer having 10 'periods'. According to Table 1 
and Fig. 2, each period consists of two sublayers - A and B. The sublayer A in each 
period maintains a constant thickness, i.e., 0.2 urn, while the thickness of sublayer B 
increases by an interval 0. 1 um from the 1st (0.2 um thick) to the 10th period (1.1 um 
thick). To account for interdiffusion induced formation of new intermetallic com- 
pounds between the coating and substrate, a transitional layer (1 um thick) having a 
gradient modulus is also introduced in the FEA simulations. 

To ensure high accuracy around the stress concentration zone, a refined mesh was 
used within an area of 20 X 20 um of the coating directly underneath the indenter 
(Fig. 7b). The total number of mesh elements is 40730, which includes the elements 
within the indenter. Further mesh refinement did not improve the simulation results 
significantly. Time-dependent deformation, such as creep, was not considered in the 
simulations. The contact between the indenter and the sample was assumed to be 
frictionless. Boundary conditions can be described as follows 31 ' 39 : the bottom of the 
simulation block is fixed in the z direction, while the right edge of the block is fixed in 
the x direction. The left edge of the block coincides with the axisymmetric axis in 
order to generate 3D results. The indenter tip was positioned at z = 0 um before 
indentation. The indentation process was simulated as the downward displacement of 
the indenter tip from 0 ~ 2.4 um at a step of 0.2 um. 

The sublayers A (made of Mo 5 Si 3 ) in the multilayer coatings are treated as an 
isotropic elastic material in simulations. In comparison, the sublayers B (dominated 
by elongated MoSi 2 grains) exhibit a columnar microstructure, and shear sliding 
would occur if the threshold stress is reached under indentation loading. 
Consequently, the sublayers B are treated as an anisotropic elastic-plastic material in 
simulations. The threshold stress governing the sliding in the sublayers B can be 
derived experimentally from the relationship of indentation hardness vs load for the 
MMC coating (Fig. 3b) 21 . It was estimated to be —9.1 GPa. Fig. 7c shows the modulus 
distribution over the thickness of the MMC coating. 
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